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Abstract 

We report on the measurement of the magnetic susceptibility and of ESR transitions in the garnet 
substance TbsGasO^ (TGG). The results are compared with a calculation in the framework of 
crystal field theory for the orthorhombic surroundings of the six inequivalent Tb ions of TGG. We 
also present a calculation of the magnetization for the three main crystal directions. 

PACS numbers: 71.70.Ch, 75.30. Sg, 76.30.-v 
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I. INTRODUCTION 



TbsGa50i2 (TGG) is a dielectric material with a cubic garnet structure. Each Tb 3+ ion 
has the same D 2 symmetry in its own local coordinate system which leads to pronounced 
crystal field (CEF) effects in various physical properties. TGG exhibits a transition to 



a antiferromagnetically ordered phase at 0.35K (ref.jl 



[2!). The terbium ions are located 



on inter-penetrating triangular lattices, but they do not form a simple Kagome lattice, 
where the triangles lie in a plane, but are arrayed in three dimensions, forming a so-called 
Hyperkagome lattice which like the Kagome lattice is geometrically frustrated. 

Various interesting experiments have been performed recently with this substance. One 
of them is the phonon Hall effect. In analogy to the classical Hall effect in conducting 
materials, the appearance of a thermal gradient in the direction perpendicular to both, the 
applied magnetic field and the thermal flux, was observed {4, 5|. Another experiment is 
the acoustic Faraday effect where the rotation of a linearly polarized ultrasonic wave in 
a magnetic field applied along the propagation direction was observed and quantitatively 



analyzed 



BE 



7| . In another experiment the temperature dependence of the symmetry elastic 



constants were measured and interpreted with a simple CEF-model 

The aim of this paper is to show, that a CEF-Hamiltonian for the orthorhombic sorround- 
ings of the Tb ions is well fitted to describe the low energy spectrum and the thermodynamic 
properties of TGG. We expect that this is also of intereset for a detailed model of transport 
phenomena like the phonon Hall effect. 

First we discuss the CEF-Hamiltonian, next we calculate the the magnetization in the 
three important crystal directions in high magnetic field and compare the result with existing 
experiments [9|. Afterwards we show experimental results of the magnetic susceptibility and 
compare it with analogous calculations. Finally we discuss our ESR results and interpret 
them using the same CEF-model. 

II. CEF-MODEL 

A Tb 3+ ion has 8 4f-electrons leading with Hunds rule to S=3, L=3 and J=6. The 
CEF for these ions with local orthorhombic symmetry can be described with a Hamiltonian, 
introduced by Guillot et al. in 1985, which in our notation reads: 
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TABLE I: Crystal field parameters bij for TGG in cm 1 as given in ref. 



7~Ll oc — B 2 qO 2 + B22O22 + -B40C4O + -£"42042 + B44O44 (1) 

+ -B60C60 + B 62 62 + B G4 Oq4 + B m O eG + gfisBJ . 



Here -B is the magnetic field in the local coordinate system of a Tb 3+ ion. The Oij are 
the Stevens operators as defined in [l(|, and the B^ are the crystal field parameters. The 
B^ can be calculated from the by given in ref. 



1,1 

9|. We list the hj from [9|] in table Hand give 



the explicit relation between the and the B^ in the appendix. For a detailed discussion 



of the conventions used in defining the crystal field parameters see e.g. [11]. 

In this paper we use the Hamiltonian given by eqfj] and take into account the six inequiv- 
alent ions in the unit cell [3] to calculate the magnetization, the susceptibility and ESR 
matrix elements of TGG. The local axis of the coordinate systems and the rotation matrices 
connecting the local and laboratory coordinate systems are given in table ILT1 

Diagonalization of the Hamiltonian eqjl] gives the energy levels for the six ions as a 
function of magnetic field. These are given in figiU For B=0 we find a low lying quasi 
doublet, 3.7K apart, followed by a quasi triplet, with energy levels between 60K and 64K 
and a singlet at 67. 8K. All other remaining singlets are above 433K. For cubic symmetry we 
would have a ground state doublet and an excited triplet followed by a singlet. The splitting 
from cubic to orthorhombic symmetry is therefore small. 

In finite magnetic fields we observe crossovers of the energy levels. As field direction we 
take the three experimentally most common cases, the magnetic field parallel to the cubic 
[110] c direction (see ref. (9)]), parallel to cubic [100] c direction (see ref. 6, 7]) and parallel to 
cubic [lll]c direction, corresponding to the main axis of the crystal. The most interesting 
crossovers occur for the lowest energy levels in figHh at 9T for B || [110] c , in figjlfe at 20T 
for B || (100) and in fig© at 27T for B \\ [lll] c . 
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TABLE II: Rotation matrices from the local coordinate systems of the six inequivalent Tb ions to 
the laboratory system (column 1). Axis of the local coordinate systems of the six ions (column 
2,3,4). 

III. MAGNETIZATION 

Since the energy levels of TGG display a complex structure of level crossings it is worth- 
while to consider the magnetization, though due to the strong magneto-caloric effect {3] , the 
magnetization of TGG in pulsed high fields is difficult to measure. We calculate 

(4) = ~Tr{J k exp(-/?(?f + gfM B Bj))} (2) 

where ( J l k ) is the contribution of the ith ion to the kth component {k = x, y, z) of 

4 




FIG. 1: Figure a,b,c show the six lowest energy levels for field in [110] c direction, figure d,e are 
the corresponding levels for field in [100] c direction, and figure f,g for field in [lll] c direction. The 
corresponding ions are marked in the figures, the energies are in Kelvin and the magnetic field is 
in Tesla. 

the magnetization in the laboratory system. The components of the total magnetization 
vector are given by (Jk) = Y^=i(^l)- ^ is the partition function. Note that the W are 
the Hamiltonians in the laboratory system and also B is now the magnetic field in the 
aboratory system. We obtained the % l by a rotation of 7-Li oc from the local systems 
lsjljljj . Equivalently one can of course calculate the magnetizations in the local systems 
and transform them back to the laboratory system, which for a simple vector quantity like 
the magnetization would be definitely the method of choice. However in view of a more 
complex study in future concerning the elastic constants of TGG in high magnetic field 
we choose to rotate the Hamiltonians. This means that the W we are dealing with in 
our calculations are more complex than the Hamiltonian given in eqJU since e.g., also 
with j odd appear. The crystal field constants in the laboratory system were calculated by 
Mathematica. Once these rotations are done the calculation of the oberservables discussed 
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FIG. 2: Ground state magnetization \(J) \ = ({J x } 2 + {J y ) 2 + (Jz) 2 ) 1 ^ 2 with the field parallel to the 
[110] c , [001] c , [lll] c . The saturation is at \(J)\ = 6. 

in this paper can all be conveniently performed in the laboratory system. 

FigfJ shows the total magnetization |(J)| = ((J x .) 2 + (J y ) 2 + (J z } 2 ) 1/2 at T=0 for the 
three field directions. The steps and plateaux in the magnetization mark cross-overs of the 
lowest levels as pointed out above. Changes of slope also occur because of the summation 
of the contributions of the different ions. So far only the magnetization in [110] c direction 



was measured and calculated. For this direction our results agree with Guillot et al. We 
want to point out, that the magnetizatons in the other directions also exhibit interesting 
cross-over phenomena, which could be used to corroborate the validity of the underlying 
CEF- model up to high fields. Crossover effects at 20T for B || [100] c have been observed 
previously in acoustic measurements 



Since in ref . {9 ] the magnetization was given in units of — ^ B , we make a quantitative 



comparison with our calculation. We get from ref. {9] for B = 15T an experimental value 
of ~ 41.5 ^ B . for 2 TGG molecules and from our calculation |(J)|[no] = 4.72/(Tb ion). 
This gives for 6 Tb ions and = 1.5 a value of 42.5 2mo ^ ulcs in close agreement with the 



experiment of ref. {9] 



6 



0.08 



0.06 
0.04 



0.02 



"0 100 200 300 

T in Kelvin 

FIG. 3: Measured magnetic susceptibility x ( re d points) and result of CEF-calculation (black line) 
versus temperature in Kelvin. The inset shows the same data plotted as 1/x compared with a 
Curie- Weiss fit (green dashed line). All the susceptibilites shown are dimensionless. 

IV. MAGNETIC SUSCEPTIBILITY 

The magnetic susceptibility was measured in the temperature range between 2K < T < 
300K and in a magnetic field of B = 0.1T using a Quantum Design SQUID magnetometer. 
All measurements were performed on a single crystal with the external field oriented parallel 
to the cubic [100] direction. The data were corrected for the temperature-independent 
diamagnetic core contribution and the magnetic contribution of the sample holder. The 
latter was determined from an independent measurement. The experimental data versus 
temperature are shown in fig|3j 

A fit by a Curie- Weiss law x^ 1 = (T - 0)/C from 10K to 300K gives 6 = -7K and 
the Curie constant C = 0.2277K. Evaluating C gives an effective magnetic moment of 
p cff = 9.25 close to p eff = J(J + 1) = 9.72 for J = 6 at high temperatures. A 6 = — 7K 
compared with a Neel-temperature T/v = 0.35K implies a strongly frustrated spin system in 
the ordered anti-ferromagnetic region. 

We also calculated the magnetic susceptibility of TGG using the CEF scheme, that is we 
calculated x l {T) of ion i in the laboratory system from the free energy with the formula 
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Density 


p = 7.22 g/cm 3 


Molecular weight 


M = 1017.37 


Molvolume 


V M = 140.91cm 3 /Mol 


Number of Tb-ions per cm 3 


n = 1.28 x 10 22 cnr 3 


Lande g-factor 


g L = 1-5 



TABLE III: Physical parameters for TGG 



X\T) = = ^ (EJn\J z \n)e W (^)) 2 

-lrfrEJH-/»| 2 exp(^) 
+ |E,, m ™ ! exp(^) with ,- = l,..6 (3) 

and obtained the total susceptibility x(T) = Ei=i X'CO by summing over the contribu- 
tions from the six inequivalent ions. It is obvious that the magnetic susceptibility of cubic 
TGG is isotropic. This can be also seen from fig. 2 where the initial slopes for the 3 differ- 
ent field directions are the same. The calculated susceptibility is plotted as a function of 
temperature together with the experimental data in figl3] . The measured Xm in cm 3 /Mol 
TGG has been divided by the molvolume using the parameters of table IIIII Without any 
fitting parameters the agreement between experiment and calculation can be considered as 
very good, except for a small deviation between 20K and 50K. 

As an example we compare the calculation with the experiment for T = 20K. From eqJH] 
we get that = 0.548T" 1 at T = 20K. With Mi = n g L /i B < J t > and the values 

of table IIIII this leads to Xm — 0.97 x 10 _2 cm 3 /cm 3 which has to be compared with the 
experimental value x cxp = 0.88 x 10~ 2 cm 3 /cm 3 . 

V. ESR EXPERIMENTS 

The ESR experiments were done employing a multiple-frequency 14 T ESR spectrometer 
[jjj] equipped with a millimeter-wave vector network analyzer (product of AB Millimetre) 
and a 16 T ESR spectrometer equipped with VDI radiation sources (product of Virginia 
Diodes Inc.) similar to that described in [if]]. The external field was applied parallel to the 
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FIG. 4: ESR transitions measured at T=2.5K, with the field parallel to the [110] c direction. The 
lines correspond to ESR modes as described in the text. 

propagation direction of the electromagnetic radiation and with perpendicular incidence on 
the [110] c face of the TGG crystal. The electromagnetic field vectors e,b are in free space 
transverse waves and the Zeeman Hamiltonian reads 

Hg = gfi B J Kf (4) 

with linear polarized waves and propagation direction /c||[110] c , which in the local coor- 
dinate systems corresponds to the z-direction and the (y|>±!>!) directions respectively. 
We performed the measurements in the frequency range 200-600 GHz. 

Experimentally we found 3 ESR lines, labelled A,B and C. FigJH shows the resonances 
we found at 2.5K and figElgives the temperature dependence of the resonance lines up to 97 
K. Trying to identify these lines within our calculated level scheme (fig. 1 a,b,c) we find that 
there are only few transitions which we can expect to be detectable and which we actually 
have observed in ESR. 

To ensure that the transitions we consider have nonvanishing weight we also calculated 
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FIG. 5: Temperature dependence of the three observed ESR-lines A,B,C at a a frequency of 209.6 
GHz. The field is applied in [110] c direction. Note, that the ESR spectra are offset for clarity. 

the ESR transition matrix elements (m\J x \n) and (m\J y \n) between levels n,m. The relevant 
transition probabilities are shown in the fig. El Since however the microwave intensities used 
in the experiment are strongly dependent on the frequency we can only make a qualitative 
comparison with the calculated matrix elements. 

In detail we considered the following transitions in the level scheme as possible candidates 
for the experimentally observed lines A,B and C: For the ions of type 1 (fig la) and magnetic 
field below 9 Tesla the transition from the quasidegenerate doublet to the second state has 
nonvanishing matrix elements (fig 6 d) (note that we count the ground state as line 0, line 1 
represents the first excited state, etc). Also transitions from exited states i.e. between line 
2 and line 3 or line 4 have nonvanishing weight (see figure 6 a and b). For ions of type 1 and 
magnetic field larger than 9 Tesla we expect that the transition between the ground state 
and line 1 or 2 should be detectable (fig. 6c), however the transition (0-2) has vanishing 
weight in this field range (fig. 6d). 
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For the ions of type 2,3,5,6 (fig lb) there are two possible ESR-transtions: For magnetic 
field below 9 Tesla the transition Ool (see fig. 6e) and the transition 1 -h- 2 (see fig.6f). 

For ions of type 4 (fig lc) we do not expect to see any transitions in ESR. 

Next comparing the calculated level scheme with the experimental results shown in fig. H] 
we come to the following conclusion: 

Line A: From the frequency and temperature dependence we can conclude that we deal 
here with a transition from the ground state. It is an increasing branch and can be observed 
in the field region between 1.7 to 4.7 Tesla (see figf4]). It can be nicely fitted within the 
levels 0-1 for ions 2,3,5,6 (figfjTJand figE]) and gives, using hv = g^sS a g-factor 9.34 from 
the level scheme and from experiment we find a g-factor of g = 9. Above 9T one cannot 
follow this line further because the corresponding energy levels (0-1) are constant and the 
matrix elements are vanishing. 

Line B: Also line B is a ground state excitation between the energy levels 0-1 and 0-2 for 
ion 1 (fig{T]a) but observed in experiment only for B > 10. 5T. Also the corresponding matrix 
element (0| J x ,y\^) is vanishing below 9T (figJH]) which thus correponds to our experimental 
finding. Note that matrix element (0|J y |2) has some weight for magnetic field close to 9T 
on both sides. Since this is a decreasing branch we could observe such a transition only 
below 7.4T with our lowest microwave frequency of 200 GHz. But for these fields the matrix 
element becomes rather small. Therefore line B can only be observed for B > 9T. The 
experiment gives data up to 14T despite the fact that the corresponding matrix elements 
become rather small. This can be due to the fact that the microwave intensity is rather 
large for the higher frequencies as mentioned above. The g factor for the observed B-line is 
again g = 9 and from our calculation we find g=8.7. 

Line C: Finally the observed line C has spectral weight above and below the level tran- 
sition of 9T as seen in figlH It has a g factor of g = 1.9 and cannot be fitted easily into 
the level scheme of figiH Since it is visible already at low temperatures we conclude that it 
cannot be a transition within excited states, in particular it cannot be identified with the 
transition 1-2 for ions of type 2,3,5,6 which has a g-factor g=1.6. 

But there are other possibilities for line C. One obvious idea is a resonance line due Tb 4+ 
impurities. It is known that the garnet structure can have Tb 4+ with L = 0,S = J = | like 
Gd 3+ . This could give a line with g = 2 as observed. By tilting the sample in the magnetic 
field by 10° and 20° respectively the line splitted up considerably for the frequency v = 96.8 
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FIG. 6: Calcluated ESR transitions matrix elements |(m|Jj|n}| 2 , Where i=x is depicted as black 
line and i=y as red dashed line. The abscissa is always the magnetic field in Tesla. 

GHz. This should not occur for a S-state ion with complete isotropy. But a splitting occurs 
for line A too. Another possibility is a Tb 3+ ion on octahedral site substituting a Ga 3+ ion 
(see ref.jsi). 

Considering the energy level diagram for the six different Tb 3+ ions and the given fre- 
quency range at our disposal the following points should be noted: We cannot observe the 
splitting of the quasi doublet with AE = 3.7K = 77GHz. We have not observed other 
transitions from excited states although we took measurements up to temperatures of 100K. 
The observed resonances are all on branches with increasing field. 

To summarize, we find that the two ESR-lines A and B shown in figJH correspond very 
well with our calculations. The effective g-factor is g = 9.34 for resonance A for B < 9T 
and ions 2,3,5,6. For resonance B and B > 9T we find an effective g-factor g = 8.7 for ion 
1. as compared to the experimental value of g = 9 for both line A and line B. 

We want to emphasize here, that the description by the model given in Eq. 1 is of course 
limited. Also the crystal-field parameters, which are usually determined at low magnetic 
fields might be slightly field dependent and it is not clear to what extend the small but 
nonzero exchange interaction modifies the energy levels. 
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TABLE IV: Coefficients /, 



VI. SUMMARY 



Taking all the presented experiments together we can say that the CEF scheme, intro- 
duced by Guillot et al {9], can explain magnetization, magnetic susceptibility and ESR lines 
quite well. It is rather satisfying that all these widely different magnetic effects presented 
in this paper originate in a complicated garnet structure with 6 Tb ions with different local 
symmetries, which can be quantitatively described. In the dicussions for the magnetization 
M, the magnetic susceptibility x an d the ESR we have neglected interaction effects between 
the ions. The existence of T/v and G point of course to interaction effects. But they are too 
small to be observed in the magnetization M and in ESR. 



VII. APPENDIX 



The Bij used in eq{T] can be calculated form the given in table U by the following 
equations: 



B 2j 
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and the are listed in table [TV] 
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